Cooperative resonance linewidth narrowing in a planar metamaterial 
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We theoretically analyze the experimental observations ol a spectral line collapse in a metamaterial 
array ol asymmetric split ring resonators [Fedotov et al, Phys. Rev. Lett. 104, 223901 (2010)]. 
We show that the ensemble ol closely-spaced resonators exhibits cooperative response, explaining 
the observed system-size dependent narrowing ol the transmission resonance linewidth. We lurther 
show that this cooperative narrowing depends sensitively on the lattice spacing. 
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Resonant multiple scattering plays an important role 
in mesoscopic wave phenomena that can also be reached 
with electromagnetic (EM) fields. In the strong scat- 
tering regime, interference of different scattering paths 
between discrete scatterers can result in, e.g., light local- 
ization [lj-an effect analogous to the Anderson localiza- 
tion of electrons in solids. Metamaterials comprise artifi- 
cially structured media of plasmonic resonators interact- 
ing with EM fields. Due to several promising phenomena, 
such as the possibility for diffraction-free lenses result- 
ing from negative refractive index there has been a 
rapidly increasing interest in fabrication and theoretical 
modeling of such systems. The discrete nature of closely- 
spaced resonators in typical metamaterial arrays raises 
the possibility to observe also strong collective radiative 
effects in these systems. 

In recent experiments Fedotov et al. observed a dra- 
matic suppression of radiation losses in a 2D planar meta- 
material array Q . The transmission spectra through the 
metamolecular sheet was found to be strongly dependent 
upon the number of interacting mctamolccules in the sys- 
tem. The transmission resonance quality factor increased 
as a function of the total number of active resonators, fi- 
nally saturating at about 700 metamolecules. The meta- 
material unit cell in the experiment was formed by an 
asymmetric split-ring (ASR) resonator, consisting of two 
circular arcs of slightly unequal lengths. The currents 
in these ASRs may be excited symmetrically (antisym- 
metrically), yielding a net oscillating electric (magnetic) 
dipole; Fig.Q] 

In this letter we theoretically analyze the collec- 
tive metamaterial response, observed experimentally in 
Ref. We find that strong interactions between a dis- 
crete set of resonators, mediated by the EM field, charac- 
terize the response of the ensemble and results in collec- 
tive resonance linewidths and frequencies. We show how 
the cooperative response of sufficiently closely-spaced res- 
onators is responsible for the observed dramatic narrow- 
ing of the transmission resonance linewidth (increasing 
quality factor) with the number of resonators |3| . In par- 
ticular, the system exhibits a collective mode with an al- 
most purely magnetic excitation, uniform phase profile, 
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FIG. 1: (Color online) Asymmetric split ring metamolecules. 
(a) The symmetric mode with the currents in the two meta- 
atoms oscillating in-phase and (b) the anti-symmetric mode 
with the current oscillating tt out-ol-phase. For the symmetric 
case the dominant contribution is a net electric dipole moment 
in the plane ol the resonator and lor the anti-symmetric case 
a net magnetic dipole moment normal to the plane of the res- 
onator, (c) The uniform phase-coherent collective magnetic 
eigenmode where all the metamolecules exhibit a magnetic 
dipole normal to the metamaterial plane, (d) Numerically cal- 
culated collective eigenmode v m (for u ~ 0.28A, 5uj ~ 0.3r). 
The amplitude and phase ol magnetic (green) and electric 
(blue) dipoles are represented by the length and direction ol 
the arrows, respectively. 



and strongly suppressed radiative properties with each 
ASR posscssessing a nearly equal magnetic dipole mo- 
ment. We show in detail how this mode can be excited by 
an incident plane wave propagating perpendicular to the 
array through an electric dipole coupling to an ASR, even 
when the magnetic dipole moments are oriented parallel 
to the propagation direction. We calculate the resonance 
linewidth of the phase-coherent collective magnetic mode 
that dramatically narrows as a function of the number 
of ASRs, providing an excellent agreement with experi- 
mental observations. The linewidth is also sensitive to 
the spacing of the unit cell resonators, with the closely- 
spaced ASRs exhibiting cooperative response, due to en- 
hanced dipole-dipole interactions. Moreover, at the reso- 
nance, and with appropriately chosen parameters, nearly 
all the excitation can be driven into this mode. Due to 
its suppressed decay rate, the transmission spectrum dis- 
plays a narrow resonance. 

Our analysis demonstrates how essential features of 
the collective effects of the experiment in Ref. Q can 
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be captured by a simple, computationally efficient model 
in which we treat each meta-atom as a discrete scat- 
tcrer, exhibiting a single mode of current oscillation and 
possessing appropriate electric and magnetic dipole mo- 
ments. Interactions with the EM field then determine 
the collective interactions within the ensemble. More- 
over, our analysis indicates the necessity of accounting 
for the strong collective response of metamaterial systems 
and interference effects in multiple scattering between the 
resonators in understanding the dynamics and design of 
novel meta-materials. Strong interactions between res- 
onators can find important applications in metamaterial 
systems, providing, e.g., precise control and manipula- 
tion of EM fields on a sub- wavelength scale jij, in de- 
velopments of a lasing spaser [sj], and disorder-related 
phenomena Q. 

In order to analyze the experimental observations of 
the transmission spectra, we provide a simple theoretical 
model that is computationally efficient in large systems, 
which are required to generate strong collective effects, 
but that captures the essential features of interactions 
between the resonators. We consider an ensemble of N 
metamaterial unit elements, metamolecules, each formed 
by n discrete meta-atoms, with the position of the meta- 
atom j denoted by r j (j = 1, . . . , n x N). This ensemble 
is driven by an external beam with electric field Ej n (r, t) 
and magnetic field B; n (r, t) with wavelengths much larger 
than the size of the meta-atoms. Consequently, we treat 
meta-atoms as radiating dipoles and ignore the higher- 
order multipole-field interactions. We assume that each 
meta-atom j supports a single eigenmodc of current os- 
cillation with the charge Qj(t) as its dynamic variable. 
Then the associated electric and magnetic dipole moment 
for the meta-atoms are dj = Qjhjdj and m.j = IjAjOij, 
respectively, where d., and nr, denote unit vectors and 
Ij(t) = dQj/dt is the current. Here hj and Aj are the 
corresponding proportionality coefficients (with the units 
of length and area) , depending on the specific geometry of 
the resonators. In the dipole approximation the polariza- 
tion and magnetization arc given in terms of the density 
of electric and magnetic dipoles P(r) = X)jPj( r ) ancl 
M(r) = J2j Mj( r )i where the polarization and the mag- 
netization of the resonator j are Pj(r, t) = dj5(r — r,-) 
andMj(r,<) = xa.j8(r— Tj), respectively. Here denotes 
the position of the meta-atom j. 

Incident EM field drives the excitation of the current 
oscillations, generating an oscillating electric and mag- 
netic dipole in each meta-atom. The resulting dipole 
radiation from the metamaterial array is the sum of 
the scattered electric and magnetic fields from all the 
meta-atoms Es(r, t) = 5ZjEs,j(r, t) and Hs(r, t) = 
SjHs,j(r,t), where Eg^r, i) and Hgj(r, t) denote the 
electric and magnetic field emitted by the meta-atom j. 
The Fourier components of the scattered fields have the 
familiar expressions of electric and magnetic fields ra- 
diated by oscillating electric and magnetic dipoles Q. 



We express these in terms of the positive frequency 
components Eg(r, t) and Hg(r, t) for the frequency £1 
(k = O/c), where the total field E(r) = E+(r) + E"(r) 
and (E+)* =E", 
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Here G denotes the radiation kernel representing the elec- 
tric (magnetic) field emitted from an electric (magnetic) 
dipole [7[ . The explicit expression for the corresponding 
radiated field from a dipole with amplitude v reads 



G(r, k) ■ v = (f x v) x f ■ 
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where f = r/r. Similarly, G x (r, fc) represents the radia- 
tion kernel for the magnetic (electric) field of an electric 
(magnetic) dipole source • Specifically, the correspond- 
ing radiated field from a dipole with amplitude v yields 

Gx(r,fc)-v=^(l- I i-)fxv. (4) 

Each meta-atom interacts with the incident fields 
Ej n (r, f) and H in (r,t), as well as the fields radiated by 
all other meta-atoms in the metamaterial sample and 
with its own self-generated radiation. In each meta-atom 
the oscillating current generates electric and magnetic 
dipoles which both radiate magnetic and electric fields, 
according to Eqs. (fTJ) and J2|). These fields couple to 
dynamical variables of charge oscillations in other meta- 
atoms, producing more dipolar radiation. In order to 
calculate the EM response of the system, we solve the 
coupled set of equations involving all the resonators and 
the fields. A system of n x N single-mode resonators 
then possesses n x N collective modes of current oscil- 
lation. Each collective mode exhibits a distinct collec- 
tive linewidth (decay rate) and resonance frequency, de- 
termined by the imaginary and real parts of the corre- 
sponding eigenvalue. The resulting dynamics resemble a 
cooperative response of atomic gases to resonant light in 
which case the EM coupling between different atoms is 
due to electric dipole radiation alone [1, [j| . The crucial 
component of the strong cooperative response of closely- 
spaced scatterers are recurrent scattering events [f| |To| 
- in which a wave is scattered more than once by the 
same dipole. Such processes cannot generally be mod- 
eled by the continuous medium electrodynamics, necessi- 
tating the meta-atoms to be treated as discrete scatter- 
ers. An approximate calculation of local field corrections 
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in a magnctodielectric medium of discrete scatterers was 
performed in Ref. ll| where the translational symmetry 



of an infinite lattice simplifies the response. 

In this work we consider asymmetric metamolecules 
consisting of two meta-atoms with different resonance 
frequencies ojj, centered around the frequency ujq, with 
\ujj — loq I <C ljq. The meta-atoms are illuminated by 
a near-resonant incident field with a narrow bandwidth 
and dominant frequency fio- For simplicity, we assume 
that the radiative electric and magnetic decay rates of 
each resonator Te and Ym are independent of the res- 
onator j and sufficiently slow Y e /m *C ojq. Then the 
dominant contribution of the dipole radiation in Eqs. Q 
and (jU is generated at the wavenumber k ~ flo/c and a 
single isolated meta-atom behaves as a simple linear LC 
circuit with ojj = 1/ y/LjCj, where we have introduced 
an effective capacitance Cj and self-inductance Lj asso- 
ciated with the circuit current mode [ijj ■ In the damping 
rate r = Ye + Ym + Yo, we add to radiative dissipation 
also a nonradiative decay Yo, resulting, e.g., from ohmic 
losses. The interaction processes between the different 
resonators, mediated by dipole radiation and originating 
from Eqs. (JTJ) and (J5J), arc analogous to frequency depen- 
dent mutual inductance and capacitance, but due to the 
radiative long-range interactions, these can substantially 
differ from the quasi-static expressions for which G x is 
also absent. 

In the transmission resonance experiment Q an ASR 
formed a unit cell resonator (metamolecule). A sin- 
gle ASR consists of two separate coplanar circular arcs 
(meta-atoms), labeled by j € and separated by 

a = r r — r;; Fig. [T] The current oscillations in each 
meta-atom produce electric dipoles with orientation d 
(d _L a) and magnetic dipoles with opposite orienta- 
tions m r = — rh; (m r _L a, d). An asymmetry between 
the rings, in this case resulting from a difference in arc 
length, manifests itself as a difference in resonance fre- 
quencies with uj r = ujq + Suj and uji = ujq — Suj. In a 
single, isolated ASR the radiative interactions between 
the two resonators results in eigenstates analogous to su- 
perradiant and subradiant states in a pair of atoms. In 
order to analyze these eigenstates, we consider the dy- 
namics of symmetric c + and antisymmetric c_ modes 
of current oscillation (Fig. [T]) that represent the exact 
cigenmodes of the ASR in the absence of asymmetry 
Suj = 0. Excitations of these modes possess respective 
net electric and magnetic dipoles and will thus be re- 
ferred to electric and magnetic dipole excitations. The 
split ring asymmetry Suj ^ 0, however, introduces an ef- 
fective coupling between these modes in a single ASR, so 
that c± = (— 7±/2 =F iA)c± — iSujc T + F±, where 7± and 
A denote the decay rates and a frequency shift, respec- 
tively, and F± represents the driving by the incident field. 
Due to the asymmetry-dependent coupling between the 
modes, the incident field with E; n || d and B; n _L rh only 
excites the electric mode when Suj = 0, but for Suj =/= it 
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FIG. 2: (Color online) Cooperative metamaterial response 
displaying collective resonance narrowing, (a) The resonance 
linewidth 7 of the collective magnetic mode v m in the units 
of an isolated meta-atom linewidth T as a function of the 
number of metamolecules N; the lattice spacings u = 1/4A 
(solid line), 3/8A (dashed line), 1/2A (dot-dashed line), and 
3/2 A (dotted line). The intermediate dot-dashed line corre- 
sponds to an asymmetry Suj = O.ir, while Suj = for all other 
curves. The lower dot-dashed line incorporates nonradiative 
loss To = 0.01T with all other curves assuming To ~ 0. (b) 
Overlap O m (b/) of v m with the state bf excited by an incident 
field resonant on the mode v m (for u = 1/2A, and Yo = 0). 
For 7 <C r there is a range of asymmetries Suj for which the 
incident field almost exclusively excites the mode v m . (c) 
Comparison between experimentally measured transmission 
resonance quality factors Q = Q/Qo (stars) from Ref. 
where Qo — 4.5 denotes the single ASR quality factor (for 
A ~ 2.7cm), and numerically calculated resonance linewidth 
7 of the collective magnetic mode v m , both with u ~ 0.28A. 
Here Suj ~ 0.3I\ Y Q ^ 0.14F 



can resonantly pump the anti-symmetric magnetic mode. 

In order to model the experimentally observed collec- 
tive response Q we study an ensemble of identical ASRs 
(with a = ae x and d = de y ; Fig. [T]) arranged in a 2D 
square lattice within a circle of radius r c , with lattice 
spacing u, and lattice vectors (ue x ,ue y ). The sample 
is illuminated by a cw plane wave Ej^(r) = ^£e y e lk ' r 
with k = ke z , coupling to the electric dipole moments of 
the ASRs. In the experimentally measured transmission 
resonance through such a sheet 0] the number of active 
ASRs, N, was controlled by decoupling the ASRs with 
r > r c from the rest of the system with approximately 
circular shaped metal masks with varying radii r c . The 
resonance quality factor increased with the total number 
of active ASRs, saturating at about N = 700. 

We find that an incident plane-wave drives all meta- 
molecules uniformly and is phase-matched to collective 
modes in which the electric and/or magnetic dipoles os- 
cillate in phase. In the absence of a split-ring asymmetry, 
only modes involving oscillating electric dipoles can be 
driven. These modes strongly emit perpendicular to the 
array (into the ±e z directions) enhancing incident wave 
reflection. The magnetic dipoles, however, dominantly 
radiate into EM field modes within the ASR plane. This 
radiation may become trapped through recurrent scat- 
tering processes in the array, representing modes with 
suppressed emission rates and reflectance, and resulting 
in a transmission resonance. In order to quantify the ef- 
fect, we study the radiation properties (for = Ym) of 
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the numerically calculated collective magnetic eigcnmodc 
v m of the system [Fig.[TJd)] which maximizes the overlap 
O m (b A ) = Wm h A\ 2 /Y,i l v f b A| 2 with the P urc magnetic 
excitation t>A = (1, -1, . . . , 1, -l) T /V2N (all ASRs os- 
cillating anti-symmetrically). We then show that the in- 
troduction of an asymmetry Suj in the resonances allows 
the excitation of v m by the incident field. This mode 
closely resembles that responsible for the experimentally 
observed transmission resonance 0, [l3[ • 

Figure^ shows dependence of the resonance linewidth 
7 of the collective mode v TO , which is the closest repre- 
sentation of the phase-coherent uniform magnetic exci- 
tation, on the number of metamolecules TV for different 
lattice spacings u (A denotes the resonance wavelength of 
v m ). In the absence of ohmic losses and for sufficiently 
small u and Suj, 7 oc 1/JV for large N. The split ring 
asymmetry only weakly affects v m . For Su) = O.OIT, the 
curves representing 7 are indistinguishable from those 
for Suj = 0. For the relatively large Suj = 0.1T, however, 
7 is increased for N > 200. The cooperative response 
and linewidth narrowing sensitively depends on u. For 
larger u (e.g., u = 3/2 A), 7 becomes insensitive to N, 
indicating the limit of independent scattering of isolated 
metamolecules and a diminished role of cooperative ef- 
fects. 

As with an isolated ASR, an asymmetry Su 7^ gen- 
erates an effective two-step coupling between collective 
electric and magnetic modes: one leg of the transition 
is provided by the coupling of the incident field to col- 
lective electric dipole excitations, while the other is pro- 
vided by the coupling introduced by the asymmetry. We 
illustrate this in Fig. [21(b) , showing the relative popula- 
tion O m (hf) of the resonantly driven collective magnetic 
mode v m , where b/ is the state induced by the uniform 
driving resonant on v m . We find that for 7 -C T and 
Suj > 7, one can excite a state in which more than 98% 
of the energy is in the target mode v m . For Suj -C 7, any 
excitation that ends up in v rn is radiated away before 
it can accumulate; the array behaves as a collection of 
radiating electric dipolcs. For larger Suj, the population 
of v m decreases since the increased strength of the cou- 
pling between the symmetric and anti-symmetric modes 
begins to excite also other modes with nearby resonance 
frequencies. Although the density of modes which may 
be excited increases linearly with N, the corresponding 
reduction of 7 means that a smaller Suj is needed to excite 
the target mode, and there is a range of asymmetries for 
which v m is excited. The narrowing in 7 combined with 
the near exclusive excitation of this mode implies that for 
larger arrays the radiation from the sheet is suppressed 
and hence the transmission enhanced as in Ref. 0. In 
Fig. EJc) we compare the experimentally observed trans- 
mission resonance [H to our numerics for To — 0.14P that 
results in the expected saturation of quality factor with 
N . We use the experimental spacing u ~ 0.28A, when 



Suj ~ 0.3r is chosen consistently with the measurements 
in Ref. (l3j . The excellent agreement of our simplified 
model that only includes dipole radiation contributions 
from each meta-atom can be understood by a notably 
weak er q uadrupole than dipolar radiation field from an 
ASR [14J . The result also confirms the importance of the 
uniform magnetic mode v m on the observed transmission 
resonance. The observed saturation is due to a combi- 
nation of a fixed Suj, which in larger arrays leads to the 
population of several other modes in addition to v m , and 
ohmic losses in the resonators which set an ultimate limit 
to the narrowing of 7. Since in the experiment Suj > 7, 
the resonance narrowing could be improved by using a 
smaller Suj, while a sufficient transfer of excitations to 
v m can still be achieved; Fig. [21(b). 

In conclusion, we analyzed the recent observations of 
transmission spectra in a metamaterial array of ASRs. 
We showed that the system can exhibit a strong coop- 
erative response in the case of sufficiently closely-spaced 
resonators. Moreover, we demonstrated how an asym- 
metry in the split rings leads to excitation of collective 
magnetic modes by a field which docs not couple directly 
to ASR magnetic moments. The excitation of this uni- 
form phase-coherent mode results in cooperative response 
exhibiting a dramatic resonance linewidth narrowing, ex- 
plaining the experimental findings Q. 
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